Absence of Long-Range Coherence in the Parametric Emission from Photonic Wires 
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We analytically investigate the spatial coherence properties of the signal emission from one- 
dimensional optical parametric oscillators. Because of the reduced dimensionality, quantum fluc- 
tuations are able to destroy the long-range phase coherence even far above threshold. The spatial 
decay of coherence is exponential and, for realistic parameters of semiconductor photonic wires in 
the strong exciton-photon coupling regime, it is predicted to occur on an experimentally accessible 
length scale. 

PACS numbers: 42.25.Kb, 42.65. Yj, 71.36.+C, 89.75.Kd 



I. INTRODUCTION 

A central concept in the theory of nonlinear dynamical 
systems is the so-called pattern formation phenomenon, 
where a spatially ordered structure appears in an other- 
wise homogeneous system when this is driven far from 
thermodynamical equilibrium^. Very important exam- 
ples of pattern formation are found in the transverse 
dynamics of nonlinear optical systems 2 ^, e.g. lasers, 
bistable devices, as well as planar optical parametric os- 
cillators (OPOs). 

In OPOs, the pattern formation phenomenon mani- 
fests itself as the appearance of two additional coherent 
beams, called the signal and the idler, originating from 
nonlinear conversion of the pump beam when its inten- 
sity is brought beyond a threshold valued. An interesting 
point of view on OPO operation is provided by the anal- 
ogy with the Bose-Einstein condensation phenomenon of 
equilibrium statistical mechanics^: as it happens for the 
matter Bose field in the condensate mode when the tem- 
perature is lowered below the condensation temperature, 
the population of the signal/idler modes becomes macro- 
scopic above the threshold and the field becomes coher- 
ent. Although many aspects of the OPO operation are 
successfully interpreted in terms of concepts of equilib- 
rium statistical mechanics^, care has to be paid as a 
non-equilibrium steady state is fundamentally different 
from a thermodynamical equilibrium state: it is in fact 
determined by a dynamical equilibrium between exter- 
nal driving and dissipation and does not follow a simple 
thermal Boltzmann lawful. 

At equilibrium, it is well known that long-wavelength 
fluctuations of both thermal and quantum origin can de- 
stroy long-range phase coherence in reduced dimensional- 
ity Bose systems, thus replacing true Bose-Einstein con- 
densates with quasi-condensate a 5 ! 10 ! 11 . First studies for 
the non-equilibrium case have been recently published 
in an optical contexliLiSii^ pointing out several analogies 
with the equilibrium case, but the peculiarities of the 
non-equilibrium case have not been fully appreciated yet 
nor experimentally investigated. 

The aim of the present paper is to obtain analyti- 
cal predictions for the long-range coherence properties of 



the parametric emission from one-dimensional photonic 
wires. Although the theory that we develop for reduced- 
dimensionality effects in non-equilibrium systems is a 
general one and can be extended to a large class of phys- 
ical systems showing non-equilibrium phase transitions 
and pattern formation, our attention will be focussed on 
a specific example of solid-state system in which the ef- 
fects we predict are experimentally accessible with the 
current semiconductor technology. 



II. THE PHYSICAL SYSTEM 

We consider planar semiconductor microcavities in the 
strong-coupling regim o 14 ! 15 , a system in which OPO op- 
eration has been recently observedi§*i£. The extremely 
large value of the excitonic nonlinearity is responsible 
for the ultra-low OPO threshold of these systems, which 
suggests an enhanced role of quantum fluctuations. By 
nanostructuring planar cavities, reduced dimensional- 
ity photonic wires can be obtained without spoiling the 
strong-coupling condition, nor the possibility of OPO op- 
eration^. All these facts make these systems an ideal 
candidate for the study of quantum effects on pattern 
formation in non-equilibrium systems, in particular low- 
dimensional ones. 

A schematic plot of a photonic wire is given in Fig^t: 
light is confined in the growth z direction and in the 
transverse y direction by respectively the DBR mirrors 
and the lateral etching of the structure. On the other 
hand, it is free to propagate along the longitudinal x 
direction. The photonic states are then classified^ by 
the longitudinal wave vector k, the subband index j and 
the polarization index a = {||, _L}, and can be used as a 
basis over which to expand the (vectorial) electric field 
operator: 
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The polarization and spatial profile of the mode 
E/cjo^Z/j z) depends on the details of the specific structure 
and has to be computed by solving Maxwell's equations 
for the specific geometry under consideration. Generally, 
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FIG. 1: (a) Sketch of a semiconductor photonic wire with 
an embedded quantum well (QW). Confinement in the z,y 
directions is provided by respectively a pair of Distributed 
Bragg Reflectors (DBR) and a lateral etching of the cav- 
ity, (b) Linear regime dispersions of signal, pump and idler 
subbands and sketch of the parametric process under ex- 
amination, (c) Exciton density in the pump and signal 
modes as a function of the pump laser intensity for a given 
hu p = 1.398725 eV. Full (dashed) lines represent linearly 
stable (unstable) solutions. Cavity parameters inspired by 
Ref|l3 fr^c = f^x = l-4eV, k z = 2uAim~ 1 , htt R = 2.5 meV, 
L y = 4/im, Kg = 5 x 10~ 6 eV/im 2 and ft7s, P ,i = 0.15 meV. 



it will show a non-trivial spatial and polarization struc- 
ture^. As a function of the longitudinal wavevector k the 
dispersion of the ja subband is approximately giveni^by 



to be constant along the x direction and to have a negli- 
gible overlap in the y direction with all subbands but for 
the p one. A possible choice to fulfill both conditions, is 
to use a wide laser beam incident on the sample along 
the yz plane at an angle chosen in such a way to match 
the quantized k y wavevector of the p subband. 

Because of the nonlinearity due to exciton-exciton in- 
teractions within the quantum well, two pump polaritons 
in the p subband are then parametrically converted into 
one signal polariton in the s subband and one idler po- 
lariton in the i subband. 



III. THE THEORETICAL MODEL 

As discussed in detail in RefQ, a simple, yet quantita- 
tive description of the dynamics of the polariton quantum 
field in a microcavity is based on a Wigner representation 
of the quantum fields in terms of classical, yet stochastic 
C-number polariton fields. As the parametric dynamics 
is concentrated in the three s,p, i subbands, it is use- 
ful to project the polariton field onto these modes, so 
to restrict our analysis to three one-dimensional polari- 
tonic fields <fi s ,p.i(x)- Within the Wigner framework^, 
the time-evolution of the fields <j) St p^{x) is described by 
stochastic differential equations, whose form in Fourier 
space is: 
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where k y ^ a and k z Jcr are the quantized photon wavevec- 
tors in respectively the y and z directions. 

The strong dipole coupling between cavity photons and 
the quantum well exciton results in a Rabi frequency flu 
larger than the damping rates so that the eigenmodes of 
the system at linear regime are polaritons, i.e. coherent 
superpositions of cavity- photonic and excitonic modes. 
Compared to the photonic one J2Jl, the dispersion of the 
exciton is negligible, so that the excitonic response can be 
safely considered as being local. As a consequence, each 
photonic subband only mixes with the excitonic state of 
wavefunction Ejy CT (y, z = zqw), zqw being the position 
of the quantum well^i. 

An example of the polaritonic dispersion Sj(k) of the 
three lowest s,p, i subbands is plotted in Fig^ for re- 
alistic experimental parameters. For each subband, only 
a single polarization state is considered here: as exper- 
imentally shown in Ref. IHI the splitting between the 
a =||, J- polarization states in photonic wires is in fact 
large enough for the parametric process to select a sin- 
gle polarization state, while the other one can be safely 
neglected. 

Polaritons are injected into the cavity by means of a 
coherent and monochromatic pump laser incident on the 
cavity at a frequency u>l close to resonance with the cen- 
tral p subband. In order to exclusively inject k = po- 
laritons into the p subband, the laser field amplitude has 
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dWj(k,t). (3) 



As discussed in the previous section, the incident laser co- 
herently drives the pump subband p with an amplitude 
Fi,. The damping rates "fj(k) accounts for radiative and 
non-radiative losses of both the excitonic and photonic 
components of the polaritons. In the region of interest, 
they have a weak dependence on k. In agreement with 
the fluctuation-dissipation theoremSi, damping is inti- 
mately connected to quantum fluctuations. These show 
up in the Wigner formalism as Gaussian, complex- valued, 
white noise terms dWj(k,t) satisfying 



dWj{k,t) dWj>(k',t) 
dWj{k,t) dW*,(k',t) 



(4) 
47rdtS(k- k')5jj>. (5) 



Expectation values of symmetrically-ordered observables 
are then obtained as averages of the corresponding C- 
number quantity over the noise dWj(x,t). 

The nonlinear coupling responsible for the parametric 
process comes from the collisional exciton-exciton inter- 
actions in the quantum well. All other nonlinear pro- 
cesses are in fact much weaker and do not contribute 
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in a significant way to the dynamics of the system. At 
the simplest level, exciton-exciton collisional interactions 
can be modelled by a polarization-independent, repulsive 
contact potential with strength g > Oi^i^. 

Under the hypothesis that the polaritonic dynamics is 
restricted to the s,p,i subbands of the lower polariton, 
the mean-field polariton-polariton interaction energy can 
be written as 

£ n ,(k,k/) = f^<S' S €(p - k)MP ~ k% (6) 

where the overlap integrals of the excitonic component of 
the transverse wave functions of the different subbands 
are defined as: 

The transverse excitonic wavefunctions <frx±(y) are P r °~ 
portional to Ej~j a {y, z = zqw) and normalized to the 
excitonic Hopfield weight of the subband: 

Jdy\^ x± (y)\ 2 = \U x \ 2 . (8) 

For the sake of notational simplicity, we have not ex- 
plicitely indicated the fc-dependence of the transverse 
wavefunction 4> x ± {y) '■ given the small range of wavevec- 
tors k involved in the parametric emission (of the order 
of the inverse coherence length l/£ c ), it can be safely 
neglected in the following. 

The generalization of (JJJ to the case of spin-dependent 
exciton-exciton interactions is straightforward: in this 
case, the scalar coupling constant g has to the replaced by 
a four-indices tenso»2ii2£ which is to be contracted with 
the four (vectorial) excitonic wavefunctions 4> X j_(y). 



IV. MEAN-FIELD 

The mean-field solution of the motion equations © 
once the noise terms are neglected has the simple plane 
wave structure: 

fi(x) = Jp~° e i{k ^-^ t+e ^ . (9) 

At low pump intensities, a stable solution with a finite 
density in the pump mode only (i.e. p° = p° = 0) ex- 
ists. The wavevector k p = 0, as well as the frequency 
uj p = ojl are fixed by the external laser; in particu- 
lar, we have chosen ojl to be just below the bottom of 
the p subband so as to avoid complicated pump bista- 
bility effects22i24iS. For increasing pump powers, this 
solution eventually becomes dynamically unstable above 
some threshold, and is replaced by a new solution with fi- 
nite signal and idler intensities p° i . The pump frequency 
and wavevector remain the same lo p = lol and k p = 0, 



while many values are possible for the signal/idler fre- 
quencies oj Sy i and wavevectors k St i, with the only con- 
straint that u>i = 2lul — lj s and fej = — k s and that the 
solution is dynamically stable. This latter condition de- 
fines a bandi of stable k s values and for each of them, the 
densities p° i and the frequencies u> Sl i can be simply ob- 
tained by substituting the ansatz I0 into the determin- 
istic part of the motion equation ©. While the phase 
6p of the pump field is fixed by the one of the incident 
laser, the signal and idler ones 9 St i remain free, with only 
a constraint on their sum 9 S + 0j. The phase rotation 
9 S ^ — ► 9 St i ± A9 is in fact a symmetry of the problem, 
a 17(1) symmetry which is spontaneously broken 26 above 
the parametric threshold by the solution ©. 

The results for the signal, pump, and idler densities 
p° _ i are plotted in Fig^ as a function of the pump 
laser intensity II = \Fl\ 2 for the specific case of a signal 
wavevector k s = 0, which is the first to become unstable 
for the chosen parameters. The results would however 
be qualitatively similar if a finite k s ^ was taken. The 
signal/idler densities p° si continuously grow from zero 
starting from the lower threshold, get to their maximum 
value, and then go back to zero at the upper threshold, 
where we are back to a pump-only mean-field solution 



V. BOGOLIUBOV THEORY 

The simplest way to include the effect of noise terms 
is to linearize the Wigner equation J3J around the mean- 
field solution © and treat the noise as a perturbation. 
This is best done using the same ansatz: 

(t>j{x,t) = y/pfi + 5 P] (x, t) e <(*j*-wi*+fli+Mi(x.*)}. (io) 

as in the generalized density-phase Bogoliubov theory 
developed for the study of quasi-condensates at equilib- 
riumiSiii. This approach does not require the existence 
of long-distance coherence, but only that density fluc- 
tuations are small \Spj\ -C p°, and the signal and idler 
phases are well defined and locally locked by the para- 
metric process \2S9 p — 59 s — 66i\ <§; 1. 

Simple Ito manipulations of @, lead to the linearized 
motion equation for the fluctuations in frequency space: 

LuU(k,u) = C(k)U{k,uj) + ^jY{k)/A W(k,oj). (11) 

The fluctuation 6-vector is defined as U(k,ui) — 
[u( + )(fc, to), u^\k, w)] T ,with the 3- vector: 

u£\k, u ) = ± i^5e m {k,u). (12) 

Z V Pra 

The index m takes here the values 1,2,3 corre- 
sponding respectively to the subbands s,p,i. The 
complex- valued noise 6-vector is analogously defined as 
W(k,Lu) = [w(fc, uj), — w*(— fc, — w)] T , and is such that 
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(wj n (fc,w)w„(fc / , w')) = 8ir 2 6 mn 6(k - k')S(uj - u') and 
(w m (fc, cj)w n (fc', a/)) = 0. The diagonal 6x6 damping 
matrix T(fc) has Tu(k) = 7;(fc; + k) for 1 < / < 3 and 
Fu(k) — 7;_3(fc;_3 — k) for 4 < I < 6, while all other 
elements vanish. The 6x6 matrix C(k) has the typical 
Bogoliubov structure 



C(k) 



( M{k) Q 
\ -Q* -M*{-k) 



(13) 



in terms of the 3x3 matrices M(k) and Q whose elements 
(for 1 < to, n < 3) are given by 

M mn (k) = [e m (k m +k)- uj m - i-i m {k m + k)/2] S m , n 

T z 9lD °m,n+r~te V PrPt 
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FIG. 2: Imaginary (a) and real (b) parts of the Bogoliubov 
excitation spectrum around the mean field steady-state. Solid 
line indicates the Goldstone branch. Pump-intensity as at the 
dotted line in Fig^j. Same cavity parameters as in Fig0 



The eigenvalues of the C{k) matrix give the excita- 
tion spectrum A„(fc) of our system around the mean-held 
steady-state solution. An example of such a spectrum 
is plotted in Fig[3 Dynamical stability of the mean- 
held solution requires that for all branches Im[Aj/(fe)] < 
for all values of k. As the mean-held solution @ spon- 
taneously breaks the U(l) signal/idler phase symmetry 
d s s — * &s,i ± of the motion equation a Goldstone 
branch is necessarily present. For clarity, it will be de- 
noted in the following by the index v = G. Its main 
property is that its dispersion Xc{k) at k = is exactly 
Xo(k = 0) = 1 . The behavior of this eigenvalue Xo(k) in 
the neighborhood of k = is crucial for the long-distance 
properties of our system. Dynamical stability and ana- 
lyticity 1 ' 2 - arguments guarantee that for k — > one can 
expand 



Im[A G (fc)] ^ -fc 2 /2a, 



(14) 



with a strictly positive coefficient a > 0: differently 
from zero-sound in equilibrium Bose gases, the Goldstone 
mode is here a diffusive, non-propagating mod e - 26 ' 28 . 

Note that this result, being based on symmetry argu- 
ments, does not depend on the restriction of our model to 



a single polarization state per subband. Explicite inclu- 
sion of all polarization states would simply increase the 
dimension of C to 12 x 12. More Bogoliubov branches 
would then appear in fig|21 but no significant changes 
would affect the nature nor the dispersion relation of the 
Goldstone mode, which (as we shall see in the next sec- 
tion) is the only responsible for the long-distance coher- 
ence properties of the emission. 



VI. FIRST ORDER CORRELATION FUNCTION 

The hrst-order spatial coherence of the parametric 
emission is determined by the spatial coherence of the in- 
cavity polaritonsL 2 ^ (added ref . to baas) . In the Wigner 
representation, this can be written for x± ^ x 2 as: 



gi 1] (xi - Xl) = ((f)*(xi,t)(j) s (x2,t)) 



w 



(15) 



Neglecting in (|1(J|> the density fluctuations Spj which do 
not contribute to long-distance properties, and expand- 
ing higher order correlations by means of the Wick the- 
orem for Gaussian stochastic variables^, one obtains: 

gi'Hx! - x 2 ) = p °e ik °^- x ^e- x *( Xl - X2 \ (16) 

where the phase-phase correlation function is defined as: 



1 



Xs(xi -x 2 ) = -{[S9 s (x u t) - Se s (x 2 ,t)] 2 ) 



w ' 



(17) 



This can be evaluated using 111(1 and ((12JI : 

(i-e*-)r„ ( ,A-, 



Xs(x) 



1 



32tt 2 p% ^ 



uj - X* u (k)][uj - X^k)Y 
(18) 



in terms of the matrix elements 

f^(fc) = P^B^k)Y{k)Bl{k)P. 



(19) 



The vector P is defined as P = (1, 0, 0, -1, 0, 0) T and 
the matrices B v (k) are the projectors on the eigenspace 
corresponding to the eigenvalue X u (k). The frequency 
integral in H18[) can then be performed closing the inte- 
gration contour in the complex plane: 



Xs{x) 



16irp 



F V / dk(\-e lkx ) — 



T utl {k) 



(k)-X,(k) 



(20) 



As usual, the large distance properties are determined 
by the k ss region, and the dominant contribution 
comes from the Goldstone mode, that is for v = p = G . 
A simple analytical expression can be obtained by keep- 
ing only this mode, then approximating the smooth func- 
tion TcG(k) with its k — value Tq G (of the order of the 
7's), and finally extending the integral to infinity. This 
results in an asymptotic exponential decrease 



(21) 
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with a coherence length l c equal to 



«r° GG 



(22) 



In FigEl £ c is plotted as a function of the pump power. 
Its behaviour closely follows the one of the signal density 
p° so that £ c goes to at the edges of the OPO region, 
and has its maximum in the middle. This prediction 
can be compared to the numerical results of a full Monte 
Carlo simulation of the Wigner equation as reported in 
Ref.7. The agreement is good and the small quantitative 
discrepancy is mostly due to the three-mode expansion 
performed here, which does not take into account the 
upper polariton. 
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FIG. 3: Analytical prediction (solid) and QMC data 

(circles) for the coherence length as a function of the pump 
power. Same cavity parameters as in Fig0 



For a given value of the mean-field interaction energy 
gpj and all detuning and cavity parameters kept fixed, 
the coherence length has the remarkable scaling £ c oc 1 / g. 
The extremely high value of the excitonic nonlinearity 
implies that the effect of quantum fluctuations is strongly 
enhanced as compared to OPOs based on conventional 
nonlinear optical medic*Sii2i2Si2£. In particular, one can 
observe in Fig|3]that the predicted coherence length is of 
the order of hundred microns for realistic experimental 
parameters, which means that it should be within experi- 
mental reach with the current semiconductor technology. 

As a final point, it is very instructive to compare the 
non-equilibrium formula ll'-l) with the one for a degener- 
ate ID Bose gas of density n at equilibrium. The sym- 
metry which is spontaneously broken is in fact the same 
U(l) in the two cases2&. While at T — the equilibrium 
system shows a power law decay of correlations, an ex- 
ponential decay is recovered at finite T, with a coherence 
length given byii: 



C q = 



2nh 2 
ksTm 



(23) 



The dependence on the density is the same in the equi- 
librium and non-equilibrium cases, while the role of the 



temperature T is played in the non-equilibrium case by 
the noise associated to damping. Another, more strik- 
ing, difference concerns the free boson mass m, which is 
replaced in the non-equilibrium case by the parameter a 
defined in (|14f> . i.e. the inverse of the second derivative 
of the imaginary part of the Goldstone mode dispersion. 



VII. CONCLUSIONS 

In conclusion, we have developed an analytical theory 
which is able to describe the long-distance phase coher- 
ence of the light emitted by a one-dimensional optical 
parametric oscillator above threshold. Because of the re- 
duced dimensionality, quantum fluctuations turn out to 
be able to destroy the long-range order. At long dis- 
tances, the coherence function shows an exponential de- 
cay with a coherence length £ c inversely proportional to 
the nonlinear coupling constant. 

This result is very general and holds for any kind of 
one-dimensional parametrical oscillator, as well as for 
other examples of non-equilibrium phase transitions in 
reduced dimensionality. We have concentrated our at- 
tention on semiconductor photonic wires in the strong 
coupling regime simply because such systems appear as 
very promising candidates to experimentally investigate 
this physics. The extremely large value of the exciton- 
exciton interactions corresponds in fact to a value of £ c 
short enough to be within reach of the current semicon- 
ductor technology. 

In addition to their interest for the fundamental 
physics of pattern formation and non-equilibrium phase 
transitions, the conclusions of the present paper may also 
have significant consequences for device applications of 
OPOs, as quantum fluctuations provide an intrinsic lim- 
itation to the long distance phase coherence that can be 
obtained even in the absence of other decoherence and 
disorder effects. 
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